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Abstract. Internucleon interactions evolved via flow equations yield soft po- 
tentials that lead to rapid variational convergence in few-body systems. 

The Similarity Renormalization Group (SRG) [2] provides a compelling 
method for evolving internucleon forces to softer forms by decoupling low- from 
high-momentum matrix elements [21 H] . A series of unitary transformations pa- 
rameterized by s (or A = s~^/^) is implemented through a flow equation: 



where T^ei is the relative kinetic energy. Applications to nuclear physics to date 
in a partial- wave momentum basis have used Gs = Tj-d [3] , so the flow equation 
for each matrix element is (with = {k\T,.ei\k) = h?k'^/m) 



-{k\V,\k') = -{ek - ek'f{k\Vs\k') + ^(efc + e^, - 2eg){k\Vs\q){q\V,\k') . (2) 



The flow of off-diagonal matrix elements is dominated by the first term, which 
drives them rapidly to zero. This partially diagonalizes the momentum-space po- 
tential, leading to decoupling [3]. Pictures showing different initial NN potentials 
evolving to band-diagonal form can be viewed at the SRG website [5j. 

In the left panel of Fig. [H the ^Sq phase shift for the Argonne vis NN poten- 
tial is shown up to 800 MeV lab energy. The phase shifts for the SRG potential 
Vs are indistinguishable at any A because the evolution is exactly unitary at the 
two-body level. To test decoupling, the original and evolved SRG potential (to 
A = 2fm~^) are smoothly set to zero for momenta above /cmax = 2.2 fm~^. The 
SRG phases are unchanged up to the corresponding -Eiabi so high momenta are 
not needed. The AV18 phases are completely changed because even low-energy 
observables have contributions from high momentum, which has led to the mis- 
conception that high-energy phase shifts are important for nuclear structure [4j. 
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Figure 1. Left: Decoupling in the ^So phase shift for the Argonne uig NN potential [4]. Right: 
Decoupling in the triton with the N^LO chiral EFT potential of Entem and Machleidt [6]. 



A similar story for the triton ground-state energy with a chiral EFT N^LO 
potential is seen in the right panel, which shows the energy as a function of 
^max- Because of decoupling, the full answer for smaller A is reached for lower 
femax- A consequence is faster convergence in variational and ab initio few-body 
calculations, as shown in the left panel of Fig. [2] for two N^LO potentials, where 
the energy is plotted against the size of a harmonic oscillator basis. A complete 
study with NN potentials in the no-core shell model is given in ref. [6]. 
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Figure 2. Left: Convergence in the triton [1]. Right: Tjon line traced out by SRG-evolved NN 
potentials, labeled by A [6]. 

The commutators in Eq. ([1]) imply that the evolving Hamiltonian will have 
many-body interactions to all orders (i.e., insert second-quantized operators). 
Thus there will always be a truncation and the evolution will only be approxi- 
mately unitary. The present calculations evolve only the NN part, which explains 
why different converged triton energies are seen in Fig. [2j This is a controlled 
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approximation in the range of A for which the variation is comparable to the 
truncation error inherent in the initial EFT Hamiltonian. The variation is seen 
to be natural in Fig. [2] and the left panel of Fig. [3l which also shows the improved 
convergence (decreasing error bars) for smaller A [6]. 
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Figure 3. Left; Ground-state energy in ^Li vs. A |6]. Right: Running of the bare (unevolved) 
deuteron quadrupole moment. 

Including the 3N interaction is essential for nuclear structure. The SRG evo- 
lution only modifies the short-distance part of the potential or operators. This 
is illustrated by the weak running of the bare quadrupole moment in the right 
panel of Fig. [3l Since the chiral EFT 3N force will be modified only at short dis- 
tance, a good first approximation should be to simply re-fit its two short-distance 
coefficients at each A. In parallel, we are implementing the 3N evolution by ap- 
plying Eq. ([T]) in the three-particle space, which does not require solving the full 
three-nucleon problem [3] . The evolution of the NN interaction is independent of 
spectators and the equation for the 3N interaction has no disconnected pieces. 
A model calculation that introduces a diagrammatic treatment is described in 
ref. [7]. 
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